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Abstract 

A long-standing puzzle in charm physics is the large difference between the D^ — )• K'^K~ 

and D^ — t- vr+vr^ decay rates. Recently, the LHCb and CDF collaborations reported a sur- 

^^ prisingly large difference between the direct CP asymmetries, /S.Acp, in these two modes. We 

^^ show that the two puzzles are naturally related in the Standard Model via s- and d-quark "pen- 

guin contractions". Their sum gives rise to lS.Acp^ while their difference contributes to the 

^ two branching ratios with opposite sign. Assuming nominal SU(3) breaking, a U-spin fit to the 

,J^ D^ — )• i^+7r^,7r+A'^ ,t:^tt^ ,K^K^ decay rates yields large penguin contractions that naturally 

^ explain AAfjp. Expectations for the individual CP asymmetries are also discussed. 
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I. INTRODUCTION 

There are several surprising experimental facts in D" decays to pairs of charged pseu- 
doscalars. The first one is the long-standing puzzle of the large rate difference between 
the singly Cabibbo-suppressed (SCS) decays, -D° — )► K^K^ and D^ — )► tt+tt^. The two 
amplitudes would be equal in the U-spin symmetric limit, whereas the measured rates yield 

A{D^ -^ K+K' 



(0.82 ±0.02), (1) 



for their CP averaged magnitudes, after accounting for phase space. ^ This has led to spec- 
ulation that U-spin breaking in SCS D decays is 0{1) [IHH], rather than of the nominal size 
characterized by 

eu ~ ifK/U - 1) ~ 0(0.2) . (2) 

However, such a conclusion is premature, as indicated by the following interesting experi- 
mental observations: 

1. The CP averaged magnitudes for the Cabibbo-favored (CF) D^ — )■ K~7t~^ and doubly 
Cabibbo-suppressed (DCS) D^ — )■ K^n^ amplitudes satisfy 



1 = (15.1 ±2.8)%. (3) 



V,dVus A{D^ -^ K-7r+) 

2. The CP averaged amplitudes satisfy the experimental "sum-rule" relation (a similar 
sum rule has been discussed in [H]) 

\A{D'> -^ K+K-)/V,sVus\ ± \A{D^ -^ tt+tt-) / ^.K^ _^ ^ (4 0±1 6 

d\ 

(4) 



— uie \j^^D.^K^^-),v^,V^^\ + \A{D^^K'^+)/V^,V^,\ ^ (4-0±1.6)%. 



The expressions in ([3]) and Q would vanish in the U-spin limit. Thus, the fact that they are 
small experimentally suggests that U-spin is a good approximate symmetry in these decays. 
The alternative is that U-spin breaking is 0{1) in the SCS decays under consideration, and 
nominal in the CF/DCS ones. However, such a hierarchy of U-spin breakings would be left 
unexplained. 



1 We define \A{D^ -> PP)\ = [T{D^ -^ PP) 87rm|,/pc]^/V(lkeV), where T is the CP averaged decay rate 
and Pc is the center-of-mass momentum of the finai state mesons. 



Another interesting result is the surprisingly large time-integrated CP asymmetry differ- 
ence, AAcp = AcpiD^ -^ K^K') - Acp{D^ -> tt+tt"), recently measured by the LHCb 
and CDF collaborations [TUl ITT] . Inclusion of the Babar and Belle measurements of the 
individual K^K'^ and tt^tt^ time-integrated CP asymmetries |12l [13] and the indirect CP 
asymmetry Ay [El [15] yields the world average for the direct CP asymmetry difference pT] 



A^^'^ = ^cp(^ ^ K^K-) - A'^JpiD -^ TT+TT-) = (-0.67 ± 0.16) %. (5) 

In the Standard Model (SM), the ratio of penguin-to-tree amplitudes is naively of 
0{}ycbVub/VcsVus\oislT^) ~ 10~^, yielding A^^p < 0.1%. This expectation is based on 
estimates of the "short-distance" penguins with 6-quarks in the loops. While AAjJjp could 
be a signal of new physics [HI [ISH23], it was argued long ago that long-distance effects could 
conceivably give large direct CP violation (CPV) due to hadronic enhancement of penguin 
amplitudes [21]. Small CPV was subsequently predicted in [25] using a model for final-state 
interactions. Making use of lessons learned from the D — t- PP data that has since become 
available, it was estimated in [26] (see also [27]) that formally power-suppressed long-distance 
s- and (i-quark "penguin contractions" can yield penguin-to-tree ratios of C(0.1%), thus po- 
tentially explaining the observed A^^p. More recent works argue that a SM explanation is 
either marginal [2B], or not possible j25] . 

In this paper we show that the possibility of a large penguin amplitude in the SM acquires 
further support from the experimental data. A consistent picture emerges in which 

1. U-spin breaking is nominal. This helps explaining (|3| and (|4]) with minimal tuning of 
strong phases. 

2. The U-spin invariant sum of the s- and c?-quark penguin contractions enhances the 
penguin amplitude, thus explaining A^^p in (|5]). 

3. The difference between the s- and c?-quark penguin contractions, which we refer to as 
the "broken penguin" (with respect to U-spin), explains the difference in the decay 
amplitudes in ([I]). 

The last point above requires a broken penguin amplitude that is of the same order as 
the tree amplitude (with a ~ 50% smaller value preferred), yielding substantial interference 
between the two (see also [7]). In turn, nominal U-spin breaking implies that the penguin 



amplitude is enhanced relative to the tree amplitude by C(l/e[/) (with ~ O.S/ec; preferred). 
This is in the favored range to explain the observed AA^p. 

The situation we describe in this paper resembles the one which arises in kaon decays: 
the apparently large isospin breaking in K — )■ tttt [30] results from a combination of nominal 
isospin breaking and the "AJ = 1/2 rule", i.e., the enhancement of the Aq amplitude relative 
to ^2. Here we suggest that the apparently large U-spin breaking in D^ — )■ K^K^ ,tt^tt~ 
decays is a consequence of both nominal U-spin breaking and an enhancement of the AU = 
penguin matrix elements relative to the tree amplitude. 

The paper is organized as follows. In Section [IT] we discuss the relevant AC = 1 effective 
Hamiltonian, explain our counting in eu for the operator matrix elements, and give the 
U-spin decomposition of the decay amplitudes. Note that our counting is modified in order 
to take explicit account of the penguin contractions, thus allowing for their enhancement. 
Implications of the measured CP averaged decay rates for the penguin contractions are then 



studied. In Section III we incorporate the CP violation data. We conclude in Section |IV[ A 
derivation of the U-spin decomposition together with U-spin breaking is given in Appendix [X| 
Appendix [B] contains the equivalent "diagrammatic" decomposition of the decay amplitudes. 
We also provide a translation between the diagrammatic and the U-spin reduced matrix 
elements. 

II. ENHANCED PENGUINS AND DECAY RATES 

The decays 

D^^K+n-, D^^K-n+, D^^n+n-, D^ ^ K+K~ (6) 

are related through U-spin. We begin with a discussion of their U-spin decomposition. Using 
Cabibbo-Kobayashi-Maskawa (CKM) unitarity we can write the Hamiltonian governing SCS 
decays as 



^SCS _ ^F 



'eff 



(7) 



-V,f,V 



cbV-ub 



1=1,2 i=3 



-Fh.c. 



where 



QPP' = {pu)v-A{cp)v-A, Qf' = iPaUi3)v-AiCfsPa)v-A (8) 



are the "tree operators", Q3,..,6 are the QCD penguin operators, and Qsg is the chromomag- 
netic dipole operator. The effective Hamiltonian for Cabibbo-favored (CF) decays contains 
only the tree operators, 

rCF _ Gf 1 






^K.K:,^agf + h.c., (9) 



and similarly for doubly Cabibbo suppressed (DCS) decays with the replacement s ^ d. 

We first work in the limit in which U-spin is exact or nearly so, but with the d and 
s quarks still of two distinguishable flavors. Our working assumption is that the penguin 
contraction contributions coming from the Q12 operators are enhanced compared to the tree 
amplitude, defined below. The enhancement is parametrized by 1/e', where e' ^ 1. Thus, 
at a scale /i ~ mo we have for the 0{l/e') U-spin invariant matrix elements 



p = {K+K-\ Y, c.Qf + Yl cmD') = 

1=1,2 j=3,..,89 

= (tt+tt-I Y C^Qf + Y C^Q.\D^) ~ 0{l/e') 

i=l,2 i=3,..,8g 

T + P = {K^K-\ Y C^QT + Y. C^Q^\D'^) = 
4=1,2 i=3,..,8g 

4 = 1,2 4 = 3,. .,83 



(10) 



(11) 



The P amplitude is a pure AU = transition. Note that the Qi,2 operators in (10) can only 
produce the K^K~ final state if the dd quark pairs are contracted (and similarly for ss in 
the case of ix^-n^). The contractions are illustrated diagrammatically in Fig. ^of Appendix 
[Bj The contractions also have a non-perturbative field theoretic definition which employs the 
lattice as a UV regulator. This is what we mean when we refer to the penguin contractions 



below. In (11) the non-contracted contributions of (5i.2 are part of T. 



The contributions of the penguin operators Q^, ..,&■, Qsg in (10) and (11) are expected to 
be an order of magnitude smaller than required to explain AA^p, see e.g., [26], and are thus 
ignored throughout this work. Note, however, that the scheme and renormalization-scale 
dependence in their Wilson coefficients cancels the scheme and scale dependence appearing 



in the penguin contraction matrix elements of (^1,2 in (10) and (11). This cancelation is 
understood whenever we refer to P below. 



The T matrix element is 0{1) in the e' counting. In the U-spin hmit it is given by 

i=l,2 i=l,2 

= (ir+TT-l Y. QQtm = {n+K-\ Y, QQf\D) ~ 0{1) . 



1=1,2 j=l,2 



This follows from the fact that both {K+tt-\, {{K+K-\ - {^T+^T-\)/^/2, {K-tt+\ and Qf , 
(<9f - Qf)/V2, Qf are U-spin triplets. 



The sum of the two amplitudes in the first line of ( 12 ) is U-spin breaking, giving a "broken 
penguin" , 

Pbreak = I ({K^K-\ Y C.(g? - Qnm + (vr+TT-l J^ C^iQf ' Qf )|5°)) . (13) 

^ 1=1,2 j=l,2 / 

The leading contribution to Pbrcak measures the difference between final-state interactions 
involving the ss and dd contractions. The broken penguin is parametrically of the size 

Pbreak ~ Cf/P ~ 0{eu/e) ~ 0{l). (14) 

In the last equality we have used the scaling e' ~ e{/, which is satisfied by the data, as shown 
below. For now, however, we keep e' and eu separate. 

As already stressed, our working assumption is that the matrix elements containing pen- 
guin contractions of (5i,2 operators are enhanced. We then have two sets of matrix elements, 
the ones that are 0{l/e') enhanced, and the ones that are not. For each of them there is also 
an expansion in the U-spin breaking parameter eu- At n-th order in U-spin breaking the 
reduced amplitudes that are not enhanced are of 0{e^), while the reduced amplitudes that 
contain penguin contractions are of 0{e^/e'). For example, summarizing the above results, 
we have the following scalings 

T~0(1), P~0(l/e'), Pbreak ~ 0(et//e'). (15) 



The expressions (12) and (13) are valid to 0(1, e(//e'). At 0{eu, e^/e') the sum of matrix 
elements in (13) also receives a contribution due to U-spin breaking in T, changing the 
l.h.s. from Pbreak — ^ -Pbreak(l " ^^sd) + ^^TiT . In principle, the two contributions - the 
-Pbrcak ~ ^uP term and the euT term - could be separated, if necessary. They correspond to 
two different topological amplitudes, with the qq fermion fields in the QW operators either 



contracted or not. In particular, the matrix elements defined in this way could be calculated 
on the lattice in the (probably not so near) future [HT] . 

The decay amplitudes are derived in Appendices [A| and [B| At order C(e(/, e^/e'), and 
using the notations of Appendix [B} they read 

AiD' -^ K^n-) =V,sV:,Til - ^6^2), (16) 

A{D' ^ TT+TT-) = - I (K.K; - V^dV:,) (r(l + ^eri) + Pbrcak(l - h^J)) 

- V:,V., (P(l - \ep) + \T) , 

A{& ^ K+K-) =\ (V;.C - KdK*d) [m - \tTi) - Pbreak(l + \ef})) 

(18j 

- K;K. (P(1 + iep) + \T) , 

A(D° ^ TT+i^-) =V,dV:,T{l + leT2). (19) 

The coefficients multiplying the U-spin breaking parameters e^ are chosen such that typically 
ej ~ 0{eu)- For simplicity the V*i^Vub suppressed terms are only given to order C(l,e(//e'), 
i.e. to first subleading order in the expansion. The scaling of the different terms in the 
V*^uh suppressed amplitudes is 

P~C(l/e'), Pep ~ C(ec//e'), T ~ C(l). (20) 

The "tree" parts of the amplitude are also given to the first subleading order in the expansion, 
which in this case is C(e(7,e^/e'), with 

T~0(1), Pbreak ~ 0(et//e'), TeTi,TeT2~0(et/), Pbreakefj ~ 0(e?./e'). (21) 
We are now ready to check how well the implicitly assumed scaling e' ~ eu compares 



with the data. In Fig. ^we display the result of a fit of (16)-(19) to the measured Z)° — >■ 
K^K~ , vr+vr", K^-a^ branching ratios [32]. Here we can safely neglect the V^iV^fe-suppressed 
contributions. There are three real parameters which are fioated in the fit - the magnitudes 
of T and Pbreak and their relative strong phase. The U-spin breaking parameters e^i, ^ti 
and e^^ are varied in the constrained range e^ G [0, 0.4] for Fig. |lj (left), and in three ranges, 
ei e [0,0.2], ei G [0,0.3], and e^ G [0,0.4] for Fig. [I] (right). The branching ratios can be fit 
perfectly, so that the x^ global minimum is x^^j^ = 0. The fit shows that 

Pbreak ~ T , (22) 

in accord with our e' ~ eu counting. This is a direct consequence of the large difference 
between the D — )■ K^K^ and D — )• tt+tt^ decay rates. 
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FIG. 1: Constraints on Pbreak vs. T from the fit to the branching ratios. Left: the contours denote 
regions allowed at la (solid), 2a (dashed), 3a (dotted). Right: the constraints at la but with e^ 
varied in the ranges e^ G [0,0.2] (solid), e^ e [0,0.3] (dashed), e^ G [0,0.4] (dotted). 

III. ENHANCED PENGUINS AND CP VIOLATION 



The time-integrated CP asymmetry for SCS D^ decays to a final CP eigenstate / is 
defined as 



Acp{D -> /) 



r{D^f)-r{D^f) 



(23) 



It receives both direct and indirect CP violation contributions (see, for example, [H]). In 
the SM the indirect CP violation lies well below the present experimental sensitivity. We 
therefore assume that within the SM the measurement of Acp{D — )■ /) equals the direct 
CP asymmetry 

I4J2_|4J2 

2rf sin 7 sin 5/ . (24) 



Af 



\Afl_-jAfl_ 



(25) 



Here we have used the fact that in the SM the CP-conjugate decay amplitudes for CP even 
final states can be written as 

Af = A{D ^ f)=A^[l + Vfe^'^^f-^^], 

Af = A{D ^f)=A^f[l + r/e'(^/+^)] . 

A'j is the dominant amphtude that is proportional to {VcsV*^ — VcdV*^), see ([7]), and r/ is the 
relative magnitude of the subleading amplitude, which is proportional to VcbV*f^. It carries 
the weak CKM phase 7 = (67.3l3;5)° [33] and the relative strong phase 5j 



'f- 



We perform a fit to the branching ratios and CP asymmetries to determine Vf for / = 
K^K~ , TT+TT^. In the fit we use the HFAG averages for the individual time-integrated CP 
asymmetries^ (which includes the Babar [12], Belle [13], and CDF measurements [31]), 



Acp{D -^ TT+TT") = (0.22 ± 0.24 ± 0.11)% , 
Acp{D -^ K+K-) = (-0.24 ± 0.22 ± 0.10)% . 

and their difference measured at LHCb [101 



(26) 



AAcp = Acp{D -^ K^R-) - Acp{D -^ tt+tt") = (-0.82 ± 0.21 ± 0.11)% , (27) 

and CDF [11] 

A^cp = (-0.62 ±0.21 ±0.10)%. (28) 

For strong phases 6f ~ 0{1) we have 

AAcp ~ 4rj, (29) 

using the fact that sin 7 ~ 0.9, and the U-spin based expectation that Acp{D — t- K^K~) 
and Acp{D — )■ tt+tt^) have opposite signs. In order to explain the central values of AAcp 
one needs 

r/ ~ 0.2% , (30) 

or, equivalent ly, 

P/T ~ 3 , (31) 

after accounting for CKM factors. In [26j the penguin contraction contributions were esti- 
mated to yield rj ~ 0.1%, (or P/T ~ 1.6), with a factor of a few uncertainty. This motivates 



us to regard a hierarchy for P/T that is much larger than (31) as unlikely. 

Our main point is that under the assumption of nominal U-spin breaking, a broken 
penguin Pbreak which explains the difference of the D^ — > K^K^ and D^ — )■ tt+tt^ decay 
rates implies a Af/ = penguin P that naturally yields ( pO| and the observed AAcp- The 
scaling Pbreak ~ ^uP together with our fit result Pbreak ~ T /2 (see Fig. [I]) yields the estimate 



rTT+-K-,K+K- 



VrrV. 



cb "ub 



v,..x. 



p 



T ± Pur. 



break 



§^-^-0.2%, (32) 



The average makes sense in the hmit of neghgible indirect CP asymmetry. 
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FIG. 2: Constraints on 5p vs. P (left) and (^p^roak ^^- -f'break (right) following from the fit to 
branching ratios and CP asymmetries at la (solid), 2a (dashed), 3a (dotted). 



for eu ~ 0.2, consistent with (30). 



In order to exhibit this result in detail, we fit the expressions in Eqs. (16)-(19) to the 



four branching ratios and the time-integrated CP asymmetries in (26)-(28). The latter are 



identified with the corresponding direct CP asymmetries, under our assumption of negligible 
indirect CP violation in the SM. This gives us eight measurements that are fitted using five 
unconstrained real parameters: the magnitudes of T, P, Pbreak and the two relative strong 

(2) 

phases. In addition there are four U-spin breaking parameters, e^i, eT2, e^d > ^p ^^^^ enter at 
the first subleading order. They are allowed to lie in the range [0, 0.4] with arbitrary strong 
phases. The x^ global minimum is Xmin — 1-14. It is not zero because the measurements of 
AAcp, -A-cpiK^K^), and ^cp(7r+7r^) are only consistent at the ~ la level. 

In Fig. [2] (left) we show the constraints on P vs. the strong phase Sp = arg(P/T) obtained 
from the fit. As expected, small values of the penguin amplitude P require a strong phase 
close to 7r/2, whereas larger values of P allow for smaller phases. It is important to note 
that the minimum value of P required at la is ~ 5.8, or roughly a factor of 2 larger than 
Tavg = 2.83, the average value of T in our normalization, which can be read off of Fig. [1} 
We also note that significant strong phases are typical for Pbreak, as shown in Fig. |2] (right), 
where constraints on Pbreak vs. the strong phase difference Sp^^^^^^ = arg(Pbreak/^) are shown. 

Our main results are contained in Fig. p^ and Fig. h We introduce the parameter e^^ , 
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FIG. 3: Constraints on P/Tavg vs. e^^ , where Tavg = 2.83. The contours denote regions allowed 
at la (sohd), 2a (dashed), 3a (dotted). 



such that 



P — f:^'-' P 

-' break t„j J . 



.-(1) 
'sd 



(33) 



.(1) ; 



as in (B12). If our fit favored e J in the nominal range for U-spin breaking, it would support 



large penguins and a SM explanation for AAcp, as in (32). 



In Fig. [3J the results for P/T^vg vs. e^^ are shown for an extended range, P < 25, 
although it is understood that the range P < 10 is physically preferred. Indeed, at 1 a we 
find that e^^ naturally falls into the nominal range [0.1,0.3], for the physically reasonable 
range P/Tavg < 3.5, equivalent to P < 10. Note that the lower edges of P/Tavg correspond 
to the minimal values of P seen in Fig. [2j Higher values of P are compensated by smaller 
values of e^^ to yield the observed difference between the K^K~ and vr+vr" rates, and by 
smaller strong phases to yield the observed CP asymmetries. 

Fig. |4] directly addresses the question of whether we can accommodate the CP asymme- 
tries with nominal U-spin breaking. In Fig. |4] we show the values for /S.Acp for the allowed 
regions in Fig. |3| for the physically more motivated range P < 10, and for an extended range 
of P. We see that we can naturally explain the world average for /S.Acp = —0.67 ± 0.16 in 
(|5| with nominal U-spin breaking. Note in particular that, while values of P > 10 allow for 
marginally larger absolute values for AAcp-, the lower bound on P arising from the need to 
explain the difference of the branching ratios translates into a lower bound on the magnitude 
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FIG. 4: Allowed values of AAcp vs. e^^ corresponding to the regions shown in Fig. 3 for P < 10 
(left) and P < 25 (right). The contours denote regions allowed at la (solid), 2a (dashed), 2>a 
(dotted). 
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FIG. 5: Result of the fit for ^cp(vr+vr-) vs. Acp{K+K-), with P < 10 (left) and P < 25 (right). 
The contours denote regions allowed at la (solid), 2a (dashed), 3cr (dotted). 



of /S.Acp (the upper edges in Fig. |4j). The somewhat unexpectedly large measured value is 
thus naturally explained. 

In Fig. p^ we show the constraints on the individual CP asymmetries, Acp{t^^t^^) and 
Acp{K^K~), that follow from our fit. In the U-spin limit we would have Acpi'n:^'^^) = 
—Acp{K^K^)- For nominal U-spin breaking and to 0{e' /eu), we have Ptt+tt- = Pr+k-- 
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FIG. 6: Result of the fit for ^c'p(7r+7r ) vs. Acp{K^K ) without the individual CPV measure- 
ments, for P < 10. The contours denote regions allowed at la (solid), 2a (dashed), 3cr (dotted). 



Thus, we expect the asymmetries to scale like 



AcpJTT+n-] 
Acp{K+K- 



A{D^ -^ TT + TT- 



-l.8{l + Oieu)) 



(34) 



A{DO -^ K+K' 

This is seen to be true for the physically more motivated range P < 10, while for larger 
values of P much larger values of \Acp{K~^K~)\ are still allowed by the data. Although 
the fit contains the individual CP asymmetry measurements as inputs, this is a non-trivial 
result given that their la intervals are substantially larger than those returned by the fit. 
For completeness, in Fig. [6] we show the result obtained without inputting the individual CP 
asymmetry measurements. A larger hierarchy for the individual CP asymmetries becomes 
possible at la. 

Finally, let us comment on the sum rule in ^, which is fulfilled experimentally to (9(4%). 
The corresponding amplitude-level sum rule (absolute values removed) is satisfied to 0{efj). 
Taking the square roots of the branching ratios at 0{eu/e'), ^ becomes 



-'sum-rule 



n 



break 



T 



Pi 



break 



T 



+ 0{eu,eye')] = 0{A%) 



(35) 



Note that Egum-ruie is not necessarily small for P\ 



break 



T and large relative strong phase. 



For instance, |Pbreak| = l^"] yields a maximum for S 



sum-rule 



of 20%, realized at 6 p. 



90°. 



However, the sum-rule decreases rapidly for smaller Pbreak and 6p. For example, for the 
choices Pbreak/^ = 0.5 (or Pbreak ~ 1-4) and 5pb^g^^ = ±45°, which lie near the two "foci" of 



13 



the la region in Fig. [2] (right), the sum rule is reduced to ~ 6.6% at 0{eu/e'), marginally 
larger than experiment. We have checked that the degree of tuning in the description of all 
observables considered is modest, about 1 part in 3, with the dominant effect due to the 



sum rule. This is what one would expect from the 0{eu, ^-fj/e') corrections to (35), assuming 
U-spin breaking of nominal size. These corrections would necessitate substantial tuning of 
Ssum-ruie if ^u Were large, thus providing further support for nominal U-spin breaking. 

IV. CONCLUSION 

We have shown that the penguin contraction matrix elements of the standard-model op- 
erators (5i,2 can provide a consistent picture for large penguins in singly Cabibbo-suppressed 
D —7- PP decays. The U-spin violating contractions of (5i,2 explain the long-standing puzzle 
of a significantly larger branching ratio for D -^ K^K~ compared to D — )■ vr+vr". At the 
same time, the U-spin conserving contractions of Qi^2 are of the correct size to naturally 
explain the large CP asymmetry difference ^Acp = Acp{D — ?• K^K^) — Acp{D — > tt+tt^) 
measured by LHCb and CDF. A crucial observation, borne out by a detailed U-spin analysis, 
is that U-spin breaking of nominal size, (9(20%), correctly relates the magnitudes of the two 
phenomena. On this basis, we conclude that large direct CP asymmetries of order a few per 
mille are not surprising given the size of Br(D — )■ K^ K^) /'Bi{D — )■ tt+tt^). 
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Appendix A: Formal U-spin decomposition 

In this appendix we perform the U-spin decomposition of the D — )■ P^P^ decays (P = 
K,7t), including U-spin breaking. The most general expressions for the amplitudes are 
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truncated at second order in U-spin breaking, that is, there are no new hadronic matrix 
elements introduced at higher orders. The U-spin decomposition is written in a form in which 
it is clear which reduced matrix elements contain the penguin contractions. We assume that 
these matrix elements are dynamically enhanced. A translation is provided between the 
amplitudes in the U-spin decomposition and the decomposition given in Section II. 

1. General discussion 

The D^ meson is a U-spin singlet. The decay operators involve a down-type quark and a 
down-type antiquark, and thus in general they can be written as a sum of a U-spin triplet 



and a U-spin singlet. The leading operators (12) form a triplet, while the U-spin singlet 
operator is proportional to VcbV*fj. It is completely negligible as far as the CP averaged rates 
are concerned, due to the CKM suppression. In terms of tensor notation the two operators 
are given by 

\ v^dv:, -\{Vcsv:,-v,,v:,)) 

(Al) 
where the {Hi) are hadronic coefficients, and Hq is proportional to the identity matrix I2x2- 
We also split the final states into a U-spin triplet and singlet 

Mi=P^ ^ ^ 1, Mo = |(7r+7r-+ir+i^-)l2x2. 

(A2) 

The U-spin breaking is induced by the nonzero strange-quark mass term, rrisss. Subtracting 
the singlet piece, the breaking introduces a spurion that is a U-spin triplet 

fe/2 \ 
M.= \ ' , (A3) 

\0 -e/2j 

where e is a small parameter which parametrizes U-spin breaking, i.e., e ~ eu- 
In the U-spin limit there are two reduced matrix elements 

tocx(/i|iJi|0), po^{fo\Ho\0), (A4) 

where /o and /i are the singlet and triplet states corresponding to Mq and Mi, respectively, 
and |0) is the U-spin singlet D^ meson. At 0{e) there are three additional reduced matrix 
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elements 

s, oc {foUH, X l,)o|0), ti oc (/i|(iJi X l,)i|0), pi (X ((A x l,)o|i/o|0), (A5) 

where 1^ represents the U-spin breaking spurion M^. At O(e^) there are three more reduced 
matrix elements, 

t2 cx ((/i X l,)o|(ifi X l,)o|0) , 4 oc ((/i X l,)i|(ifi X l,)l|0) , 

(A6) 

P2 oc (/o|i/o X (1, X l,)o|0). 

In terms of the tensor notation, we have the identities 

(iJi X l,)o = {Hu MJ = e {H,)\{V,,V:, - KdKd) h.2 

Thus, the anticommutator (commutator) of Hi and M^ projects onto a singlet (triplet) 
operator. Similarly, we have 

(/i X l,)o ={M„Mi} = e\{K^K- -Ti+T,-) h^2 

I Ti+K-\ (A8) 

(/iXl,)i=[M„Mi]=e 

l-7r-fs:+ / 

and the anticommutator (commutator) of Mi and M^ projects onto a singlet (triplet) state. 
We can now write the two-body decay Hamiltonian to O(e^) as 

Hi =toTr(/7iMi) + itiTr ([/7i,M,]Mi) + ^t^Tr ({i7i,Mj{M„ MJ) 



+ |t'2Tr ([iJi, M:^[M,, Ml]) + siTr ({i^i, MJMq), 



(A9) 



for the U-spin triplet operator, and 



n^ =PoTt{HoMo) + ^piTt{Ho{M,,Mi})+P2Tt{HoM^Mo). (AlO) 

for the U-spin singlet operator. 

Choosing a convenient final state phase convention, the decay amplitudes can be read off 



from (|A9j), (|A10j), yielding 

AiD' ^ K+TT') = V,,V:, (to - Ihe + \t'^e^) , 

A{D' ^ TT+TT-) = -\{V,sV:, - V^dV:,) (to + sie + l^e^) - V,,V:, [p, - \pie + \p2e^) , 
A(D° ^ ir+ir-) = \{V,sV:, - K,lCd) (^0 - sie + ^t^e^) - K^C (Po + ^i^ + W^) ' 
A(D° ^ TT+i^-) = KdKTs (^0 + ^tie + W^) , 

(All) 
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where we have made the replacements ti{Hi) — )■ tj and PiiHo) — )■ pi for the reduced matrix 
elements. 

The U-spin expansion is just a basis rotation, so that in general there can only be six 
independent decay amplitudes, four associated with decays to the triplet final state /i and 
two with decays to the singlet final state /q. Thus, of the three reduced matrix elements 
introduced at O(e^), only one combination is a new linearly independent amplitude. In fact, 
we can see directly from the above amplitude expressions that p2 can be absorbed into po, 
and that one linear combination of ^2 and t2 can be absorbed into to 5 leaving another linear 
combination of the two as the new linearly independent amplitude. 

There are a number of conclusions that one can draw from the above decomposition at 
different orders in e. In the U-spin symmetric limit we have the following known relations: 
(i) all 4 decay rates are equal (up to CKM prefactors); and (ii) the direct CP asymmetries 
in the SCS decays are equal in magnitude and opposite in sign. 

Working to 0{e), and neglecting the small terms proportional to VcbV*fj, the four ampli- 
tudes depend on three reduced matrix elements. Thus, there is one relation among the 1)° 
decay amplitudes, which is given by 

AK-n+ Ar+tt- _ Ak+K- ^tt+tt- ( ^^0'\ 

and similarly for the CP conjugate decays. This sum rule is broken at O(e^), whereas the 
individual amplitudes are 0{1). The experimental relation Q is also satisfied to first order, 
that is, Ssum-mie = C(e^). 

The rate difference between the K^tt^ and K^tt^ modes and the rate difference between 
the K^K~ and tt+tt" modes arise at order 0{e). However, the latter is observed to be 0(1), 
while the former is small. The immediate conclusion is that 

sie > tie. (A13) 

There are two ways in which this relation could be realized. The first one is that there is a 
U-spin breaking hierarchy which remains unexplained, i.e., a hierarchy of e's, such that the 
breaking is much larger for si than for ti. The second possibility, which is the one we have 
pursued in this paper, is that all U-spin breaking is of nominal size, but that Si is enhanced 
relative to ti due to the penguin contractions. 
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2. The penguin contractions 



Our working assumption is that at a given order in e the contractions of the ss and dd 
fields give the dominant effects. The reduced matrix elements to which the contractions 
contribute are the ones in which both the transition Hamiltonian and the final state can be 
written as U-spin singlets. Thus, they are identified with the following traces 



Tr [(Mo or {M,,M,}) x {Ho or {H,,M,})] , 



(A14) 



According to (A9) and (AlO), these are the Hq matrix elements pQ,pi,p2 and the Hi matrix 



elements Si,t2. We elaborate below, and check the consequences for the U-spin decomposi- 
tion. 

At 0(e°), Hq gives rise to the reduced matrix element Pq. This involves matrix elements 
of Q^ + Qft for the singlet final state {\K+K-) + \7t+7t-))/V2, see (|AlJ and Q. Therefore, 



Po contains both contracted and non-contracted contributions. Also at 0{e^), the f/3 = 
component of Hi gives rise to the reduced matrix element to- It involves matrix elements 
of (55 - Qf% for the triplet final state {\K+K-) - \TT+n-))/V2. Therefore, the ss and dd 
contractions must cancel in to at C(e°). 

To account for the dominance of the contractions we introduce a second small parameter 
e', such that any matrix element to which they contribute is enhanced by 0{l/e'). Without 
loss of generality, we can define it as 



Po 



(A15) 



Thus, in po the ratio of non-contracted to contracted contributions is also of 0{e'). 

The contractions enter the matrix elements of Hi for the first time at 0{e). This happens 
for the Si matrix element for which the final state is a singlet. This is because the operator 
{Hix1^)q = {Hi,M^} contains the sum Qfg+Qii) see (A7) and (13). Thus, the contractions 



do not cancel, but add up. This is also the case for the contractions in pi at 0{e), and the 
contractions in p2 and t2 at O(e^). Note, however, that pi and t2 involve the f/3 = triplet 
final state. We can transform it to a singlet with the aid of the U-spin breaking spurion i.e., 
{M^,Mq}. We are led to the following e' counting rule: a reduced matrix element that is 



associated with one of the traces in (A14) is 0{l/e'), and is C(l) otherwise. Exphcitly, 



si ~ t2 ~ Po,i,2 ~ 0{l/e'), to,i ~ t'a ~ 0(1). 



(A16) 
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(A17) 



Taking e' ~ e ~ e[/, the different amplitudes are thus 

Po~ £>(!/£), to ~ Sie ~pie~ 0(1), 

Keeping only terms to order (9(e) ~ (9(e^/e') ~ 0(e{/), i.e. the leading and subleading terms 
for both CKM structures, we finally have 

^(5° ^ K+n-) = V,sV:, (to - |tie) , 

A(D' ^ ^+^-) = -1{V,,V:, ~ V,,V:,) (to + s,e + it2e') - KbC (Po - ^i^) ' 
A(D° ^ i^+i^-) = i(K,V;; - V,,V:,) (to - sie + Ihe') - K^KTb (Po + ^Pie) , 
A(L)° ^ TT+K-) = KdC (^0 + Ihe) . 
Note that this is the most general decomposition, that is, all the subleading terms that we 



(A18) 



neglected can be absorbed into terms we kept. The decomposition in (A18) is equivalent to 



the decomposition given in (16)-(19), with the following translations: at 0{l/e), and also 



including the 0{1) non-contraction term in po, 



p^ = P+\T- 



at 0(1), and also including the 0{e) non-contraction term in si. 



(A19) 



Pie 



Pep, 



to — T, Sie — Pbrcak + o^TlT '■, 



(A20) 



and at (9(e) 



tie = Te 



T2, 



t2e^ 



-Tbreaktg^ . 



(A21) 



Fitting (A18) to the four branching ratios, dropping the VcbV*^^ terms, and taking e G 
[0, 0.4] yields the fit for sie vs. to shown in Fig. ItI It is similar to the fit for Pbrcak vs. T in 
Fig. [T| as one would expect expect, and confirms that Sie ~ to and e' ~ e. 



Appendix B: Diagrammatical representation 



Finally, we provide a derivation of the U-spin decomposition presented in (16)-(19) in 



terms of the diagrammatical representation. This will simplify the comparison with the 
result of Ref . [26] , and provides further insight into the origin of the scaling e' ~ ejy . 



19 




FIG. 7: The results of the fit of the U-spin decomposition in (A18) to the D -^ K+vr", K+K', 



vr+vr ,TT^K branching ratios. The contours denote the regions aUowed at la (sohd), 2a (dashed), 
3(T (dotted). 

The Z)° decay amplitudes Af are split into "tree-level" amplitudes AT and penguin 
amplitudes Af. The former are sums of a "tree" diagram T and an "exchange" diagram E, 
and are given in full generality by 

A^. = -\ (K,K; - V.,V:,) (T.. + E..) , 
\ (KsKTs - V^,V:,) {Tkk + Ekk) , 
VcdV*s{T^K + E^k)- 






A^ 



fBll 



4^ 



The T amplitudes are those with a u spectator quark, and the E amplitudes are the anni- 
hilation topology diagrams (cf. Fig. [s]), in which the initial cu quark pair annihilates into 
sd, ss, dd, and ds quark pairs, respectively. 

In SCS decays we can divide E and T into two contributions: those that do not involve 
the ss or dd penguin contractions, denoted Tf and Ef, and those which do, denoted Pj 
and PT and shown in the last two diagrams of Fig. pi The penguin contractions arise from 
Ql% ~ Qi% ^^"i ^^6 U-spin violating. At the quark level they correspond to the difference 
between the rescattering contributions of the ss and dd quark pairs for a given final state. 
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FIG. 8: From left to right: tree level exchange topology diagram, Ef, and the penguin contraction 
diagrams Pj and P^. 



We thus have 



-L KK — J^uu — 1 rj^ i^ — r\ 



KK 



[P. 



KK 

EA 



:), 



T,s- 
KK)i 

-iE.s^ 



T = T'^ -\- ( P^''^ P^'^\ 

■Ti" niT ' Xtttt ■kit J: 



Ekk - E^KK - {Pkk - Pkk) 



Ej, 



jpd _|_ ( pE,d _ pE,s\ 



(B2) 



where the superscripts s and d denote the identity of the relevant operators Ql\ o^ Q12 i'^^ 
the identity of the contracted quark pair in the penguin contractions). In the CF and DCS 
decays T and E of course do not receive contributions from the penguin contractions. 

Formally, the P^ contain leading power as well as power correction contributions, while 
the P^ are pure power corrections. The scheme-dependent coefficients and leading log// 
scale dependence entering the penguin contraction matrix elements cancels in the differences 
Pj'^ _ pj^^ and Pf'' - Pf''. 

In the penguin amplitudes the penguin contractions correspond to the sum of the rescat- 
tering contributions of the ss and dd quark pairs in {QI2 + Qi2)/'^- Including the am- 
plitudes for the penguin operators Q3,..,6, Q?,gi denoted P, and the above contracted and 
non-contracted amplitudes, we have 



4 = -Vc,V:,{Pf + [Pf' + Pj'1/2 + [Pf^' + Pf '^]/2 + [T] + E)]/2), 



(B3) 



where q = s{d) for / = K^K^ (tt+tt^). The scheme and scale dependence in the penguin 
contractions is now canceled by P/. However, since the contributions of the P/ to the direct 
CP asymmetries are much smaller than observed, they are neglected in our fits. 



To simplify the comparison with the previous section and with Section III, let us define 
the total non-contracted amphtudes. 



Tk-k = Tk-k + Ek-k, 






rjid 



Et 



TnK = Tt,K + Et,k 

Tkk = T^KK + Ekk > 



(B4) 
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and the total penguin contraction contributions in the tree and penguin amphtudes, (/ = 
KK (tttt) for the K^K^ (tt+tt^) final state) 

p* = pj'^ _ pf^ + pf^'i _ pf'^^ pp = (pj'<^ + pj'^ + pf'd + pf'^)/2 . (B5) 

In the latter it is understood that scale and scheme dependence has been subtracted out by 
Pf. We can now write the decay amplitudes in terms of the above quantities as 

AiD' ^ K+K-) = \ (KsKT. - V.,V:,) {Tkk - PIck) - K^KT. {Pkk + V^ + Tkk/2) . 

(B6) 

In the U-spin limit, the non-contracted amplitudes satisfy 

Tkk = %tt = Tkw = %K = T , (B7) 



where T is defined in ( 11 ) -( 12 ) in terms of operator matrix elements. Introducing the U-spin 



breaking parameters en and er2, we can express these amplitudes at 0{eu) as 

UK = ni-leTi), %i = Til + leTi) 
TK. = T{l-\eT2), %K = T{l + \eT2)- 
Similarly, in the U-spin limit, the penguin contractions in the penguin amplitudes satisfy 

Kk = K.^P, (B9) 



where P is defined in (10). Introducing the U-spin breaking parameter ep, we can write 
them at 0{eu) as 

P^^ = P(l + iep), P^^ = P(l - lep) . (BIO) 

While the Vj vanish in the U-spin limit, at 0{eu) we have 

KK = 'PL^P^re.^, (Bll) 



where Pbreak is defined in terms of operator matrix elements in (13). Pbreak and P can be 
related by 

Pbreak = e2P, (B12) 

22 



where the U-spin breaking parameter e^^ accounts for the difference in ss and dd rescat- 
tering contributions at 0{eu). The difference between Vj^j^ and P*^ enters formally at efj. 

(2) 

Introducing a a new U-spin breaking parameter e^^ to take this difference into account yields 
the relations 



V 



KK 



-Pbreakfl + qC 



2 sd ) ' 



V. 



t 

TTTT 






(B13) 



Finally, substituting above expressions (B8), (BIO), (B13) for 7/, "P?, Pj-, respectively, in 



(B6) and neglecting the P/, yields the diagrammatic expressions for the decay amplitudes 



given in (16)-(19) of Section III We have assumed that the penguin contractions are dy- 



namically enhanced by (9(l/e'), where e' ~ ejj which leads to the scalings for the various 



amplitude contributions given in (20) and (21). 
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